Quantum feedback control of a solid-state two-level system 
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We have studied theoretically the basic operation of a 
quantum feedback loop designed to maintain the desired 
phase of quantum coherent oscillations in a two-level system. 
Such feedback can suppress the dephasing of oscillations due 
to interaction with environment. Prospective experiments 
can be realized using metallic single-electron devices or GaAs 
technology. 



The principle of feedback control is used in a wide va- 
riety of physical and engineering problems. In particular, 
it can be applied in a straightforward way to tune the os- 
cillation phase of a harmonic oscillator in order to achieve 
a desired synchronization with some reference oscillator. 
An intriguing and fundamental question is whether con- 
tinuous feedback can be used to control quantum sys- 
tems; for instance, if it is possible or not to tune the 
phase of quantum coherent (Rabi) oscillations in a two- 
level system (TLS). 

At first sight the quantum feedback seems to be im- 
possible because according to the "orthodox" collapse 
postulated the quantum state is abruptly destroyed by 
the act of measurement. However, as was shown two 
decades ago, in particular by Leggett,cl in a typical solid- 
state setup the collapse of a TLS state should be consid- 
ered as a continuous process rather than as instantaneous 
event. The reason is typically weak coupling between the 
quantum system and the detector and also the finite noise 
of the detector, so that it takes some time until accept- 
able signal-to- noise ratio is reached and the measurement 
can be regarded as completed. 

While the Leggett's theory as well as the majority 
of similar approaches can describe only ensembles of 
quantum systems, the theory describing the gradual col- 
lapse of|-.a-, single solid-state TLS was developed only 
recentlyjju (A similar problem in optics was,solved much 
earlier - see, e.g. Refs.BQ and references inQ.) Basically, 
the theory says that the evolution of a single quantum 
system due to continuous measurement is governed by 
the information continuously acquired from the detect, 
tor. Similarly to classical probability, the Bayes formulaEl 
which naturally takes into account incomplete informa- 
tion from the detector, can still be applied to the density 
matrix of the measured quantum system; thus the for- 
malism is called Bayesian.u 

In case of a poor detector the extra noise acting on 
the input disturbs the measured system stronger than 
the limit determined by the uncertainty principle; this 
leads to gradual decoherence of the measured system. In 
contrast, when measured with a good (quantum- limited) 
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FIG. 1. Schematic of the quantum feedback loop main- 
taining the quantum oscillations in a two-level system (TLS). 



detector, the quantum system does not loose the co- 
herence (even though the quantum state evolves ran- 
domly)' moreover, its density matrix can be gradually 
purifi.edu that basically means acquiring as much infor- 
mation about the system as permitted by quantum me- 
chanics. 

Since the Bayesian formalism allows us to monitor the 
continuous evolution of a quantum system in a process 
of measurement, this naturally gives rise to a possibility 
of continuous feedback control of a quantum system. In 
this paper we wilLstudy the operation of a feedback loop 
proposed in Ref.u and designed to maintain a desired 
phase of quantum coherent oscillations in a solid-state 
TLS. (Quantum feedback in optics-ias been proposed 
and studied earlier - see, e.g., Refs.Q'ErEl.) In particular, 
we will study the dependence of the loop operation on 
the feedback factor F, the available bandwidth t" 1 , and 
the dephasing rate 7 e due to environment. 

As an example of the measurement setup (Fig. |l|) we 
consider a TLS represented by a single electron in a dou- 
ble quantum dot (DQD), the location of which is mea- 
sured by a quantum point contact (QPC) nearby in a way 
used in Refill. If the electron is in the dot 2 (state |2)) 
which is closer to QPC than dot 1, then the QPC tunnel 
barrier is higher and so the average current 1% through 
QPC is smaller than the average current I\ correspond- 
ing to electron in the dot 1 (state |1)). Consequently, 
from the QPC current one gets information about the 
electron location. We consider a realistic case of weak 
response, A7 = h - h < ^0 = {h + h)/2. In this 
case the measurement time S7/2(A/) 2 , which is neces- 
sary to achieve signal-to-noise ratio equal to 1 (here Si 
is the shot noise of the QPC current), is much larger 
than e/Io, so the QPC current I(t) is continuous on the 
measurement timescale and we do not need to consider 
individual tunneling events in QPC. 

The evolution of the TLS density matrix p during the 
measurement process is. described within the Bayesian 
formalism by equationsBu 
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where e and 77 are, respectively, the energy asymmetry 
and tunneling strength of the TLS [the TLS Hamiltonian 

is Htls = (e/2)(c{ c i - 4 C 2) + #(cl c 2 + 4 c i)]> and 1 = 
Jd + 7e is the dephasing rate due to detector nonidcality 
and coupling with environment. Theoretically, 7^ = 
when TLS is measured by a QPC; however, if instead of 
QPC we use a single-electron transistor, then dephasing 
7d is always significant .Q 

Notice that the ensemble dephasing rate T = 7 + 
(A/) 2 /4Sj is larger than 7 because of different evolution 
of the ensemble members due to random I(t). Individual 
realizations can be simulated using the formula 



J(i)-Io = (pH-/> 22 )AJ/2 + £(*), 



(3) 



where £(£) is the pure white noise with spectral density 
Sj = Si. If Eqs. (0)-(^) are averaged over £(£) (we use 
Stratonovich definition for stochastic differential equa- 
tions), then we get usual ensemble-averaged equations 
for TLS evolution (terms proportional to AI will disap- 
pear and 7 will be replaced by T). 

It is natural to characterize the effect of extra de- 
phasing 7d by the detector ideality, (efficiency) r\ = 
l/[l + ld AS I /(AI) 2 ]. One can showlffl that 77 = (h/2e d ) 2 
where td is the energy sensitivity of the detector includ- 
ing the effect of backaction but neglecting the correlation 
between output and backaction noises. So, an ideal case 
T) = 1 corresponds to a detector with quantum-limited 
sensitivity. 

To realize a feedback loop (Fig. fy), we can monitor the 
TLS evolution using the detector current I(t) plugged 
into Eqs. (Q)-@. Then the TLS state is compared with 
the desired state, and the difference signal is used to con- 
trol the TLS parameters H and/or e. In the example 
studied in this paper the feedback loop is designed to 
stabilize the quantum oscillations of the state of a sym- 
metric TLS (e = 0), so the desired evolution is pn(t) — 
l-pn(t) = [l + coB(m)]/2, p 12 (t) = p* 21 {t) = isin(fit)/2, 
where the frequency is SI = (AH 2 + e 2 ) 1/2 /h = 2H/H. As 
a difference ("error") signal we use the phase difference 
A(j> between the desired value 4>o(t) — ilt and the mon- 
itored value <j>(t) = arctan{2Impi 2 (i)/[pn(i) — P22 (£)]}■ 
This difference is used to control the TLS parameter H 
(controlling the barrier height of DQD); we assume a lin- 
ear dependence: Hjb = H[l — F x A<fi], where F is the 
dimensionless feedback factor. _ 

In this paper we neglect additional time delaya in the 
feedback network, however, we take into account the fi- 
nite bandwidth of a line carrying detector current (that is 
the critical parameter for a possible experiment). More 
specifically, we average the current I(t) with a rectan- 
gular window of duration r Q , I a (t) = t^ 1 J fr I(t')dt', 
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FIG. 2. Correlation function K z {t) of the TLS quantum 
oscillations for C = 1 and feedback factors F = (thin solid 
line), 0.05 (thick solid line), and 0.5 (dashed line). Nondecay- 
ing oscillations are due to synchronization by the feedback. 



before plugging it into Eqs. (|i])-(g), so that the "avail- 
able" density matrix p a (t) differs from the "true" density 
matrix p(t). Also, to compensate for the corresponding 
implicit time delay, we use Ac/) = (\> a — Q(t — KT a ) with 
K = 1/2 (we tried various k and found that k = 1/2 
provides the best operation of the feedback loop) . 

Let us start with the case of ideal detector, rj = 1, 
absence of extra environment, 7 e = 0, and infinite 
bandwidth, r a = 0. Figure @ shows numerically calcu- 
lated correlation function K z {t) = (z(t + r)z(t)) where 
z = pii — P22, for several feedback factors: F = 0, 
0.05, and 0.5. rffhe- curves are obtained using Monte 
Carlo simulationticl of the measurement process for mod- 
erately weak coupling between TLS and the detector: 
C = a(An 2 /S'/ J ff = 1 (notice that the Q-factor of 
oscillationsEj is equal to 8/C, so C — 1 is still a weak 
coupling). In absence of feedback (F = 0) the corre- 
lation function decays to zero (Fig. @) while for finite 
feedback factor the correlations remain for indefinitely 
long time (of course, assuming perfect reference oscillator 
which determines desired evolution). The nondecaying 
correlations show that the quantum feedback loop really 
provides the synchronization of quantum oscillations. 

The degree of synchronization depends on the feed- 
back factor F. One can see that for a moderate value 
of F = 0.5 the synchronization is already very good [the 
ideal case would be K z (t) = cos(Or)/2]. For the case of 
moderate or good synchronization (C/16F < 1, C/8 <C 1) 
we have derived the following analytical expression: 
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The correlation function Kj(t) = (I(t + r)/(i)} of the 
detector current I(t) have somewhat similar dependence, 
however, it also has the decaying contributionliij due to 
correlation K Z £ and (^-function contribution due to de- 
tector noise. The analytical result, 



Ki (t) = ^5(t) 
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FIG. 3. Synchronization degree D as a function of feedback 
factor F for several values r a of detector signal averaging: 
Ta/T = 0, 1/3, and 2/3, where T = 2?r/n. Dashed line 
D = exp(— C/32F) almost coincides with the upper curve. 
Dotted line corresponds to "naive" feedback with r a = T/10. 
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agrees well with the Monte Carlo results. 

The spectral density Sj(cj) of the detector current can 
be obtained as a Fourier transform of Kj(t). While in 
absence of feedback the quantum oscillations in TLS can 
provide only a moderate peak of Si(u>) around frequency 
f2 (the peak height cannot be larger than 4 times the 
noise pedestallij), the feedback synchronization leads to 
the appearance of a ^-function at the frequency of desired 
oscillations. (In principle the desired frequency can differ 
a little from f2; however, in this case the performance of 
the feedback loop worsens.) 

Besides the correlation function and spectral density, 
we have studied one more characteristic, D, of the syn- 
chronization degree. We define D as the average scalar 
product of the unity-length vector on the Bloch sphere 
corresponding to the desired state and the vector cor- 
responding to the actual state of TLS. The equivalent 
definition is D = 2(Tr pp^) — 1, where pd is the density 
matrix of the desired pure state. Perfect synchronization 
corresponds to D = 1. It is simple to show that in the 
case of weak coupling and unshifted desired frequency, 
D 2 /2 is equal to the nondecaying amplitude of K z (r) 
dependence. 

Upper solid line in Fig. || shows the dependence of D 
on the feedback factor F for C = 1 and r a = 0. One can 
see that D is proportional to F for small F ( "soft" onset 
of synchronization) and D is asymptotically approaching 
1 at large F. Using Eq. (0) it is simple to obtain expres- 
sion D = exp(— C /Z2F) which is very close to numerical 
results for moderate and good synchronization (dashed 
line in Fig. [|). 

Finite available bandwidth of the detector current I{t) 
(finite averaging time r a in our formalism) worsens the 
performance of the quantum feedback loop. The solid 
lines in Fig. show the dependence of the synchroniza- 
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FIG. 4. Dependence D(F) for C = 1, T a — 0, and several 
magnitudes of dephasing due to environment: d e = 0, 0.1, 
and 0.5. Dashed and dotted lines correspond to d e = and 
limitation of Hfb by and H/2, respectively. 



tion degree D(F) for T a /T = 0, 1/3, and 2/3, where 
T = 2ir/£l is the oscillation period. Obviously, a sig- 
nificant information loss occurs when r a becomes com- 
parable to T, leading to a decrease of D. The curves 
D(F) saturate at large F allowing us to introduce the 
dependence D max (T). Calculations for parameters of 
Fig. H show pretty good synchronization, D max — 0.993, 
for r a = T/30, while D max = 0.98 for r a = T/10, 
D max = 0.92 for r a = T/3, and D max = 0.57 for 
Ta = 2T/3. 

The main potential practical importance of the quan- 
tum feedback is the ability to suppress the effect of TLS 
dephasing caused by interaction with environment (see 
Fig. ||). Solid lines in Fig. |j show the dependence D(F) 
for several magnitudes of the dephasing due to environ- 
ment, d e = 0, 0.1, and 0.5, where d e = 7 e /[(A/) 2 /4S7] 
is the ratio between TLS coupling to the environment 
and to the detector (we still assume an ideal detector). 
First of all, we see that feedback still maintains the phase 
coherence of TLS for infinitely long time. However, for 
finite d e the degree of synchronization D saturates at the 
level less than unity. The dependence D ma x{d e ) is linear 
at small d e , numerically we obtained D max ~ 1 — 0.5d e 
for parameters of Fig. m 

Notice that the solid lines shown in Figs. [| and ^ are 
calculated assuming the TLS feedback control Hfi, — 
H[l — F x Ac/)] even when Hfb becomes negative (this 
is also an assumption for analytical results). To elimi- 
nate this unphysical assumption we have also performed 
numerical calculations with restriction Hfb > and with 
restriction Hfb > H/2. This leads to rather minor modi- 
fications of the presented curves (dashed and dotted lines 
in Fig. show the results for d e — and r a = 0). How- 
ever, important difference is that D(F) goes down at 
large F, so the optimum D max is achieved at some finite 
value of F. More detailed study of this problem will be 
presented elsewhere. 

Besides the discussed feedback based on Ac/) calcula- 



tion, we have also studied much simpler feedback loop in 
which H fb {t)/H - 1 = F{2[I a (t) - I ]/AI - cos [ft (i - 
r Q /2)]} sin[ft(i — T a /2)] (in an experiment it would elim- 
inate the necessity to solve fast the Bayesian equations). 
Quite surprisingly, such "naive" feedback can also pro- 
vide a good phase synchronization of quantum oscilla- 
tions if F/C is close to 1/4 (see dotted line in Fig. g). 
However, it requires more careful choice of F and r a than 
for the Bayesian feedback, and also suffers more signifi- 
cantly from the restriction on Hfb variation. The results 
in more detail will be discussed elsewhere. 

Experimentally, besides the realization of quantum 
feedback control of a DQD continuously measured by 
a QPC, one can also think about the TLS based on 
a single- Cooper-pair box measured by a single-electron 
transistor (see discussion inta). This realization can be 
preferable because of a rapid progress of metallic single- 
electronics technology. However, the problems are high 
output impedance of the single-electron transistor and its 
significant nonideality as a quantum detector. The third 
potential realization can be based on SQUIDs. For any 
realization the major problem is bandwidth: the feedback 
should be faster than the TLS dephasing due to environ- 
ment. Because of that, the quantum feedback of a TLS 
should probably be attempted only after the realization 
of recently .proposed Bell-type two-detector correlation 
experiment E-3 
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